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Low-dimensional systems are beautiful exam- 
ples of many-body quantum physics [1]. For one- 
dimensional systems [2] the Luttinger liquid ap- 
proach [3] provides insight into universal prop- 
erties. Much is known of the equilibrium state, 
both in the weakly 0, d, 0, 0] and strongly jl, [1] 
interacting regime. However, it remains a chal- 
lenge to probe the dynamics by which this equi- 
librium state is reached [10]. Here we present a 
direct experimental study of the coherence dy- 
namics in both isolated and coupled degenerate 
Id Bose gases. Dynamic splitting is used to cre- 
ate two Id systems in a phase coherent state [11]. 
The time evolution of the coherence is revealed in 
local phase shifts of the subsequently observed in- 
terference patterns. Completely isolated Id Bose 
gases are observed to exhibit a universal sub- 
exponential coherence decay in excellent agree- 
ment with recent predictions by Burkov et al. 
|l2l |. For two coupled Id Bose gases the coherence 
factor is observed to approach a non-zero equilib- 
rium value as predicted by a Bogoliubov approach 
|l3l |. This coupled-system decay to finite coher- 
ence is the matter wave equivalent of phase lock- 
ing two lasers by injection. The non-equilibrium 
dynamics of superfluids plays an important role 
in a wide range of physical systems, such as su- 
perconductors, quantum-Hall systems, superfluid 
Helium, and spin systems fl4 ■ [l5 ■ Il6| . Our exper- 
iments studying coherence dynamics show that Id 
Bose gases are ideally suited for investigating this 
class of phenomena. 

The starting point of our experiments is a Id Bose 
gas of a few thousand atoms trapped in a highly elon- 
gated, cylindrical magnetic microtrap on an atom chip 
[ITl llSj with typical transverse and longitudinal oscilla- 
tor frequencies of v±^ ^ 4.0 kHz and ~ 5 Hz. Our 
trapped Bose gas is in the Id quasi-condensate regime 
[H, characterized by both the temperature T and chem- 
ical potential /i fulfilling k^T^fi < hv^ [Tol ]. 

After the initial preparation of this single Id system, 
we perform a phase-coherent splitting along the trans- 
verse direction by means of a radio- frequency (rf ) induced 
adiabatic potential Uj . As shown in figure [U the final 
system consists of two Id quasi-condensates in a verti- 
cally orientated double- well potential . They are sep- 




FIG. 1: Schematic of the experiment. A single Id quasi- 
condensate is phase coherently split using rf potentials on an 
atom chip. A combination of two rf fields allows balanced 
splitting in the vertical direction ^20]. After the separation, 
the system is held in the double-well potential for a variable 
time t and is then released from the trap. The resulting in- 
terference pattern is imaged along the transverse direction 
of the system. Thermal phase fluctuations in the two quasi- 
condensates can be directly observed as local shifts of the 
observed fringe profiles. 



arated by a tunable potential barrier, whose height is 
controlled by the applied rf fields. 

This splitting process initializes the system in a mutu- 
ally phase coherent state. Directly after the splitting, the 
phase fluctuation patterns of the two individual quasi- 
condensates are identical, resulting in a vanishing global 
relative phase. This is a highly non-equilibrium state 
of the split system, which will relax to equilibrium over 
time. 

To study this time evolution of the phase coherence, 
the two Id Bose gases are held in the double- well config- 
uration for a varying time t before they are released and 
recombined in time-of-flight. The resulting interference 
pattern is recorded using standard absorption imaging 
along the transverse direction of the system. The spatial 
variation of the relative phase between the two quasi- 
condensates translates directly into local phase shifts in 
the interference pattern (Figure [2|) . 



2 



If the two parts of the system are completely sepa- 
rated (compare methods), the equilibrium state consists 
of two uncorrelated quasi-condensates. Consequently, we 
observe an increasing waviness of the interference pattern 
with time, which in the end leads to a complete random- 
ization of the relative phase 0{z,t) (Figured^, b). This 
change in the interference pattern is a direct visualization 
of the dynamics of the phase fluctuations. Qualitatively 
similar behavior was recently observed at MIT [lH for 
elongated condensates with fi ~ 2hi'± and T ~ 5hv±. 

For a finite tunnel coupling (compare methods) be- 
tween the two systems, we also observe an increase in 
the waviness of the interference. However, in contrast to 
the completely separated case, the final equilibrium state 
shows a non-random phase distribution (Figure [2t,d). 
This is caused by the phase randomization being coun- 
terbalanced by the coherent particle exchange between 
the two fractions. The final width of the observed phase 
spread depends on the strength of the tunnel coupling 



t = 1 ms 



t = 7 ms 



t = 10 ms 
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A quantitative measure of the fluctuations of the local 
relative phase is given by the coherence factor 



(1) 



where L is the length of the analyzed signal. 

We obtain the coherence factor from a single image 
by extracting the local relative phase 6{z, t) from the 
interference pattern in each vertical pixel slice. To ac- 
count for variations in the Id atomic density ni^{z) due 
to the longitudinal confining potential, we use only the 
central 40% of each image in our analysis. Over this 
range, nid{z) varies by only ^ 15% from its peak value. 
In the following, we neglect this modulation and assume 
a homogeneous density, obtained by averaging over the 
density profile in this center region. 

We first investigate the time evolution of the coherence 
factor "^{t) for the case of completely separated Id quasi- 
condensates, separated by a potential barrier sufficiently 
high to suppress any tunnel couphng (Figure [2K,b). The 
time evolution of the measured coherence factor for six 
different combinations of initial temperature Tin, Id den- 
sity rtid, and transverse trap frequency v^^ is displayed 
in figure [3l In a time window of 1-12 ms we observe a 
universal sub-exponential decay of the form 
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with decay time constants ranging from 5 ms to 9 ms 
(Table m. 

Our results can be directly compared with a recent the- 
oretical study by Burkov et al. |12l | based on a Luttinger 
liquid approach [3] , which predicts a = 2/3 for the expo- 
nent in Eq. [2l For all the cases studied, we find a to be 
between 0.64 and 0.67 (Table[l|, in very good agreement 
with the theoretical prediction. This agreement is even 



t = 4 ms 
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FIG. 2: Direct observation of the phase dynamics through 
interference, (a) Example images of the observed interference 
patterns for hold times t = 1,4,7,10 ms in the case of iso- 
lated Id systems (taken from the T = 133 nK data set in 
figure [Sja) . White lines indicate the bright nodes of the lo- 
cal interference pattern, extracted from fits to the profiles in 
each vertical pixel slice. The randomization of the relative 
phase 9{z, t) over time leads to an increased waviness of the 
observed fringe patterns, (b) Corresponding polar plots of the 
local relative phases from the central regions of five images. 
For the isolated systems we find a complete randomization of 
the local relative phase over time, (c) For a finite tunnel cou- 
pling through the potential barrier, the phase randomization 
is counterbalanced by particle tunnelling, resulting in a sup- 
pression of the fringe waviness (images correspond to figure 
llljc)). (d) This results in a non- random steady-state phase 
spread, whose width depends on the tunnel coupling and the 
temperature of the system. 



Tin 

[nK] 


nid 
[l/^m] 


[kHz] 


vi_ 
[kHz] 


a 


to 

[ms] 


Ti 
[nK] 


82(28) 


20(4) 


0.7(1) 


3.3 


0.64(8) 


9.0(4) 


76(10) 


133(25) 


34(5) 


1.2(2) 


3.3 


0.65(7) 


5.5(3) 


145(13) 


171(19) 


52(4) 


1.8(1) 


3.3 


0.64(4) 


6.4(3) 


186(15) 


81(31) 


22(4) 


0.9(2) 


4.0 


0.65(3) 


8.1(2) 


85(10) 


128(23) 


37(4) 


1.5(2) 


4.0 


0.66(3) 


5.9(2) 


153(13) 


175(20) 


51(5) 


2.1(2) 


4.0 


0.64(6) 


6.1(4) 


194(17) 



TABLE I: Measured exponents q, decay time constants to and 
final temperatures Tj after the splitting for the data shown 
in figure [3l 
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FIG. 3; Time evolution of the coherence factor for uncou- 
pled Id quasi-condensates. Double logarithmic plots of In(^) 
as a function of time t for transverse confinements of (a) 
i'± = 3.3 kHz and (b) i'± — 4.0 kHz, respectively. Three dif- 
ferent data sets corresponding to different initial temperatures 
and line densities are shown for each configuration. Each data 
point is the average of 15 individual measurements. Error 
bars indicate the standard error of the mean. In this repre- 
sentation, the exponent of the decay function is given directly 
by the slope of the observed linear evolution of ln(\l/). The 
lines shown are linear fits to the data, the obtained slopes are 
0.64±0.08, 0.65±0.07, and 0.64±0.04 for (a) and 0.65±0.03, 
0.66 ± 0.03, and 0.64 ± 0.06 for (b). All slopes show good 
agreement with the theoretical prediction of 2/3. 



more remarkable because the approximations used in ref. 

to obtain Eq.[2]and a = 2/3 are only valid for t > to, 
whereas our data covers the full range from up to 2 tg. 

At our typical temperatures (Tin ~ 80 — 150 nK), the 
dynamics of the coherence factor is dominated by ther- 
mal phase fluctuations [12, 24]. Burkov et al. derive 
the decay of '^{t) by assuming that the system of two 
Id Bose gases can evolve into thermal equilibrium. The 
sub-exponential coherence decay then results from the 
fact that damping in a Id liquid at finite temperature 
is always non-hydrodynamic, which in turn is caused 
by the break-down of superfluid order in Id on length 
scales longer than the temperature-dependent correla- 
tion length [2^. The experimentally observed decay 
agrees with this non-trivial theoretical prediction, which 
is strong evidence that the final equilibrium state of our 
system is truly thermal equilibrium. In turn, this sug- 
gests the non-integrability of the experimentally realized 
system [26 1. 



Following ref. [1^ the decay time constant to can be 



related to the parameters of the Id quantum gas: 
to = 2.61 ngnuK/Tl 



(3) 



where K = ttH 

weakly interacting Id Bose gas [27]. g 



2id. ig the Luttinger parameter for the 

2hv±as is the 



effective Id coupling constant with Og being the s-wave 
scattering length, m the mass of the atoms, and Tf the 
final equilibrium temperature of the split system. 

T[ can be extracted from the decay constant to follow- 
ing Eq. [3] using independently measured values of v± and 
7iid. The results, together with an independently mea- 
sured initial temperature of the unsplit system Tin from 
time-of- flight images, are compiled in table HI Within the 
error bounds, the temperatures and Tf agree. Never 
the less, our data seems to indicate an increasing differ- 
ence between Tin and Tf for increasing initial tempera- 
ture. 

This observation is in qualitative agreement with ref 
[l^ . which predicts the heating during the splitting to 
scale with \/TuJ(. This is also supported by the observa- 
tion that for identical Tn, we find the same temperature 
increase in the split system for both of the two different 
transverse confinements used in the experiment. 

One should note here that the temperature Tf is the 
temperature of the longitudinal excitations and is mea- 
sured on a time scale much shorter than the characteristic 
time given by the sound propagation through the sam- 
ple, or a thermalization timescale, which should be even 
longer. To check the universality of our experiments, we 
cut each sample in two and analyze each half separately. 
We obtain the same results within the statistical uncer- 
tainties. 

We now turn to the case of (weakly) coupled Id quasi- 
condensates realized when the potential barrier between 
the two Id systems is only slightly larger than the chemi- 
cal potential fi, and tunnelling between the two split sys- 
tems is significant. In this case we observe a qualitative 
change in the behavior of the coherence factor ^(i), as il- 
lustrated in figure m After an initial decay, the coherence 
factor saturates well above the value for random phases. 
This observation is a clear manifestation of a balancing 
between phase locking due to coherent tunnelling and the 
local phase fiuctuations in the Id quasi-condensate [l3| . 
It is the matter wave equivalent of injection locking two 
lasers. 

For a quantitative analysis, we use the results of ref. 
[l3| for the energy spectrum and the mode functions and 
express the final equilibrium coherence factor as 



= vfg X exp - 
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where 1 in our case) is the contribution of the quan- 

tum fluctuations, and 
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FIG. 4; Time evolution of the coherence factor for cou- 
pled Id quasi-condensates. (a) After an initial decay we ob- 
serve a time-independent coherence factor determined by the 
strength of the tunnel coupling. Different colors correspond to 
the different transverse double- well potentials shown in (b-d) 
as indicated. Each data point is the average of 12 individ- 
ual measurements. Error bars indicate s.e.m. The light grey 
points are the T — 128 nK data set from figure |3]d, showing 
the sub-exponential decay in the case of uncoupled conden- 
sates for comparison. The corresponding transverse poten- 
tial is shown in (e). The dotted black line in (a) indicates 
the Rayleigh test confidence level above which the considered 
phase distribution of a single measurement can be identified 
as non-random with a 90% probability. In the coupled case, 
the final phase distribution in each realization can be clearly 
considered as non-random. 



oscillation frequency [1; 



h ■ 



(6) 



^From our measured 7 we estimate « 27r x 80, 200, 900 
Hz for the three presented data sets, which should be 
compared to the coherence decay time scale of the order 
of 5 ms (Figured]). This suggests that the balancing pro- 
cess between coherent tunnelling and phase fluctuations 
stabilizes on timescales of the order of a single or few 
Josephson oscillations in a coupled Id system. 

In summary, we have studied the non-equilibrium 
phase dynamics in a Id Bose gas with the aid of mat- 
ter wave interferometry. For completely separated sys- 
tems we observe the sub-exponential coherence decay 
predicted by Luttinger hquid theory For coupled Id 
Bose gases, the balance between decoherence and phase 
stabilization due to tunnelling establishes a finite coher- 



ence 
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I . In both examples our experiments illustrate 
how the thermal equilibrium state of a phase-fluctuating 
Id quasi-condensate [ES,SS01 is approached from the 
highly non-equilibrium situation of two phase coherent 
copies created by splitting a degenerate Id Bose gas. 

The above investigation demonstrates that with the 
exquisite control provided by atom chips, one can now 
study in detail the non-equilibrium dynamics of coher- 
ence in a single quantum degenerate Id system. This 
allows us to address the fundamental issue of thermaliza- 
tion in (almost) integrable systems [i^. Our approach 
can easily be extended to study the dynamics of two- and 
three-dimensional quantum systems, and to include other 
factors such as internal state dynamics. Moreover, it will 
be the starting point for detailed studies of coupled Id 
quantum systems. One prominent example is the quan- 
tum Sine-Gordon model, which plays an important role 
in many fundamental physics questions [3, 0, 13, lH . 



are the energy eigenvalues and amplitudes of the Bogoli- 
ubov modes, respectively. 7 thereby quantifies the tunnel 
coupling between the two Id systems. Independent mea- 
surements of T and riid allow us to extract 7 directly 
from the measured equilibrium coherence factors. For 
the red, blue, and black data sets in figure IH we obtain 
■^^j-^ ~ 0.003,0.01, and 0.2, respectively. Comparison 
with numerical calculations of the tunnel coup ling , based 
on the spatial overlap of the wave functions [28], yields 
qualitative agreement within a factor of three. This sug- 
gests that a two-mode model for the tunnelling dynamics 
is not sufficient to describe our experiments. Alterna- 
tively, if 7 can be determined independently, the equi- 
librium coherence factor can be used for thermometry in 
the coupled system 



23|. 



Particle exchange between the two condensates takes 
place on a timescale given by the inverse of the Josephson 
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Methods 

Preparing a Id BEC on an atom chip 

We start by preparing a magnetically trapped ultra 
cold [T ~ sample of Rb^^ atoms in the \F ~ 



2, mp = 2 > state on the atom chip [17|, [18| using our 



standard procedure of laser cooling, magnetic trapping 
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and evaporative cooling [34I. The sample is then trans- 
ferred to a highly elongated magnetic trap {i/± ~ 4.0 kHz, 
longitudinal confinement i/^ < 5 Hz) at a distance of 
75 /im from the atom chip surface and cooled to quan- 
tum degeneracy by tuning the evaporative cooling radio 
frequency (rf) in between the ground and first excited 
transverse trapping state. The resulting sample an ef- 
fectively one-dimensional 0, with /i, fceT < hv^ con- 
taining 1...10 X 10^ atoms. 



Smoothness of the trapping potentials 

Our atom chip wires [s^ provide exceptionally smooth 
trapping potentials [3J|. The residual roughness due to 
current-flow perturbations in the trapping wire is much 
too small to be measured directly at our operation dis- 
tance h — 75 /im from the atom chip surface. An esti- 
mate based on the current flow pattern in the wire recon- 
structed from magnetic field microscopy measurements 
employing Id condensates at short distance (~ 10 ^m) 
[3a | shows that for length scales below 40 /im, the remain- 
ing potential perturbations are smaller than 1 (corre- 
sponding to < lO^"^/!). Due to the exponential damping 
of the higher spatial modes (wave vector fc = ^) of the 
magnetic field variations with distance h from the chip 
surface following exp{~kh), the disorder potentials will 
be even smaller on the short length scales corresponding 
to the phase coherence length or the healing length in 
the Id Bose gas. 



Preparing two phase-coherent Id systems 

To prepare two mutually phase coherent Id systems, 
we employ a combination of static and rf magnetic fields 
on our atom chip. Coupling electronic ground states of 
the magnetically trapped atoms results in dressed-state 
adiabatic potentials, whose versatility stems from the de- 
pendence of the potential on the angular orientation be- 
tween the rf field and the static trapping field [sH]. In 
particular, one can create very robust double-well po- 
tentials, which allow the coherent splitting of a trapped 
Bose condensate jl^ . In the experiments presented here, 
we employ a setup similar to that in ref. [2^ where the 
combination of two rf fields generated by auxiliary wires 
on the atom chip allows the realization of a compen- 
sated double- well potential in the vertical plane. For this 
double-well orientation, the observed interference fringes 
in the atomic density are horizontal, parallel to the atom 
chip surface. This allows us to image the interference pat- 
tern along the transverse direction of the system (Figure 



Preparing two Id systems with a variable tunnel 
coupling 

To introduce a variable tunnel coupling between the 
two Id systems, we adjust the potential barrier between 
3 and 8 kHz by changing the amplitude of the rf fields 
(Figure m. The actual barrier height is determined from 
spectroscopic measurements of the rf potentials to a pre- 
cision < 1 kHz 37[ . To evaluate whether the two samples 
in the resulting double-well potential can be considered 
as one-dimensional, we compare the chemical potential 
H and the thermal energy ksT of the trapped ensembles 
to the single-particle level spacing in the double-well. In 
the case of large potential barriers, the transverse con- 
finement of each individual well can be approximated by 
a harmonic oscillator, and the level spacing is given by 
the oscillator frequency h'±. For small barriers, this ap- 
proximation fails and one has to numerically calculate the 
single particle states in the transverse double- well poten- 
tial. For the configurations shown in figure |4] (b-d), we 
find energy separations between the ground-state doublet 
and first excited states of A = 2.8,3.4,3.8 kHz, respec- 
tively. Consequently, /x, /cbT < hA for all configurations, 
justifying the Id treatment of the individual systems. 



Measuring the interference pattern and extracting 
the local relative phase 

We observe the interference pattern created by the two 
expanding, overlapping atomic clouds using standard ab- 
sorption imaging. Employing diffraction-limited optics, 
optimized light path and pulse duration, and a weak mag- 
netic field to define a quantization axis, we achieve atom 
shot noise limited imaging with 3 fim resolution and a 
noise floor of ~ 3 atoms per 3x3 fim pixel. We extract 
the local relative phase 9{z) by fitting a cosine function 
with a Gaussian envelope to the observed density distri- 
bution in each individual vertical pixel slice. The free 
parameters of these fits are the relative phase 9, the con- 
trast, and the fringe spacing. The width, amplitude, and 
center position of the envelope are determined indepen- 
dently from a Gaussian fit to the total integrated density 
pattern of the central area of each image. From the mea- 
sured phases, the coherence factor is evaluated (compare 
eq. [p. This has to be contrasted to the methods used in 
ref. |23l . [s^ , where the interference patterns are summed 
up and the contrast is analyzed. 



Evaluating the coherence decay constant 

Several aspects must be considered when evaluating 
the decoherence constant to. First, the time t — 
when the tunnel coupling vanishes and the two conden- 
sates start to evolve independently, changes for different 
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atomic densities and trap parameters and has to be evalu- 
ated for each data set separately. Higher atomic densities 
[3^ , or stronger transverse confinement lead to increased 
tunnel coupling as the overlap of the two matter waves in 
the tunnelling region becomes larger [40| . The moment 
of the decouphng of the two systems can be estimated by 
numerically calculating the time-dependent tunnel cou- 
pling throughout the splitting. Experimentally, one ob- 
serves a qualitative change in the detected interference 
patterns [HI, 0, EH when the barrier becomes higher 
than the chemical potential. In the case of the widely 
split independent condensates, the uncertainty in the de- 
termination of t = is much smaller than 1 ms. 

The second aspect to consider is that for a finite sys- 
tem length L, the coherence factor does not approach zero 
but converges to a finite value. The average coherence 
factor calculated from a limited number of samples is 
non-zero, even if these phases are totally random. Such 
a limited number of samples for a system of length L 
is given by the finite imaging resolution, and by the fi- 
nite phase-coherence length itself 4] , which in our case is 
comparable to a single pixel width. This results in an ad- 
ditional offset in the equilibrium coherence factor, which 
can be determined from the number of data points used 
in the calculation of 431. 
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